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Abstract 

The nonextensitivity parameter q occuring in some of the applications of 
Tsallis statistics (known also as index of the corresponding Levy distribution) 
is shown to be given, in q > 1 case, entirely by the fluctuations of the param- 
eters of the usual exponential distribution. 
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There is an enormous variety of physical phenomena described most economi- 
cally (by introducing only one new parameter q) and adequately by the so called 
nonextensive statistics introduced some time ago by Tsallis [JJ . They include all sit- 
uations characterized by long-range interactions, long-range microscopic memories 
and space-time (and phase-space as well) (multi) fractal structure of the process (cf. 
[|TJ for details). The high energy physics applications of nonextensive statistics are 
quite recent, but already numerous and still growing, cf. Refs. @, 0, |[ ^ [| |7], || . All 
examples mentioned above have one thing in common: the central formula employed 
is the following power-like distribution: 
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which is just a one parameter generalization of the Boltzmann-Gibbs exponential 
formula to which it converges for q — > 1: 



G q =i = g(x) = c ■ exp 



x 
A 
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When G q (x) is used as probability distribution (Levy distribution) of the variable 
x G (0, oo) (which will be the case we are interested in here), the parameter q 
is limited to 1 < q < 2. For q < 1, the distribution G q (x) is defined only for 
x G [0, A/(l — q)]. For q > 1 the upper limit comes from the normalization condition 
(to unity) for G q (x) and from the requirement of the positivity of the resulting nor- 
malisation constant C q . However, if one demands in addition that the mean value 
of G q (x) is well defined, i.e., that (x) = A/(3 — 2q) < oo for x G (0, oo), then q is 
further limited to 1 < q < 1.5 only. In spite of numerous applications of the Levy 
distribution G q (x), the interpretation of the parameter q is still an open issue. In 
this work we shall demonstrate, on the basis of our previous application of the Levy 
distribution to cosmic rays @, that this Levy distribution G q (x) (P emerges in a 
natural way from the fluctuations of the parameter 1/A of the original exponential 
distribution (Q) and that the parameters of its distribution /(1/A) define parameter 
q in unique way. 



Let us first briefly summarise the result of |J. Analysing experimental distri- 
butions dN(x)/dx of depths x of interactions of hadrons from cosmic ray cascades 
in the emulsion chambers, we have shown that the so called long flying component 
(manifesting itself in aparently unnexpected nonexponential behaviour of dN(x) / dx) 
is just a manifestation of the Levy distribution G q (x) with q = 1.3. This result must 
be confronted with our earlier analysis of the same phenomenon [[| . We have demon- 
strated there that distributions dN(x) / dx can be also described by the fluctuation of 
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u = >rj 1 > 0-2 (3) 



the corresponding cross-section a = A j- (where A denotes mass number of the 
target, is the mass of the nucleon and A is the corresponding mean free path). 
The fluctuation of this cross-section (i.e., in effect, fluctuations of the quantity 1/A) 
with relative variance 

<^ 2 >_^> 2 

allow us to describe the non-exponentiality of the experimental data as well as the 
distribution G q= i,s(x) mentioned above. We therefore argue that these two numer- 
ical examples show that fluctuations of the parameter 1/A in the g(x; A) result in 
the Levy distributions G q (x; A). 

Actually the above quoted example from cosmic ray physics is not the only one 
known at present in the field of high energy collisions. It turns out J3], ||] that 
distributions of transverse momenta dN{pn r )/dpT are best described by a slightly 
non-exponential distribution G q {pt) of the Levy type with q = 1.01 -j- 1.2 depend- 
ing on situation considered. The usual exponential distribution dN(pT)/dpT = 



g{p T ) ~ exp(— ym 2 + p^/kT) contains as a main parameter the inverse tempera- 
ture (3 = 1/kT and the above mentioned numerical results leading to G^i.oih-i^Pt) 
can again be understood as a result of a fluctuation of inverse temperature j3 in 
the usual exponential formula g(pr)- This point is of special interest because of 
recent discussions on the dynamical possibility of temperature fluctuations in some 



collisions, cf. Ref. (T^, [11], |12| . Later on we shall use it to illustrate our results 
concerning q. 

To recapitulate: we claim that (for q > 1) the parameter q is nothing but a mea- 
sure of fluctuations present in Levy distributions G q (x) describing particular pro- 
cesses under consideration. To make our statement more quantitative, let us analyse 
the influence of fluctuations of parameter 1/A which are present in the exponential 
formula g(x) ~ exp(— x/X) on the final result. Our aim will be a deduction of the 
form of the function /(1/A) which leads from an exponential distribution g(x) to 
power-like Levy distribution G q (x) and which describes fluctuation about the mean 
value 1/A , i.e., such that 

G.fc Ao) =C.(l + ±±y = C, f exp (-£) / (I) A (i) (4) 

where for simplicity we have introduced the abbreviation a = From the repre- 



sentation of the Euler gamma function we have [13 



i + f = ^ / '/cc 1 < • 
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Changing variables under the integral in such a way that = j one immediately 
obtains eq. (U) with /(1/A) given by the following gamma distribution 




and variation 




Notice that, with increasing a variance decreases and asymptotically (for a — > oo, 
i.e, for q — > 1) the gamma distribution (0) becomes a delta function 8(X — A ). The 
relative variance (cf. eq.(^)) for this distribution is given by 




We see therefore that, indeed, the parameter q in the Levy distribution G q (x) de- 
scribes the relative variance of the parameter 1/A present in the exponential distri- 
bution g(x; A). 

Some remarks on the numerical results quoted before |2|, || are in order here. 
Notice that the value of q = 1.3 for cosmic ray distribution dN(x)/dx obtained 
in p| leads to the relative variance of the cross section uj = 0.3 whereas in we 
have reported value uj' = 0.2. This discrepancy has its origin in the fact that in 
numerical calculations in || we have used a symmetric Gaussian distribution to 
decribe fluctuations of the cross section, whereas the relation (||) has been obtained 
for fluctuations described by gamma distribution. In the gaussian approximation 
we expect that 

< J < q-1, (10) 

q 2 

where lower and upper limits are obtained by normalizing the variance of the /(1/A) 
distribution to the modial (equal to (2 — q)/Xo) and mean (equal to 1/A ) values, 
respectively. Therefore for q = 1.3 one should expect that 0.18 < uj' < 0.3, which is 
exactly the case. 

Let us now proceed to the above mentioned analysis of transverse momentum 
distributions in heavy ion collisions, dN(p T )/dpT 0. It is interesting to notice that 
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the relatively small value q ~ 1.015 of the nonextensive parameter obtained there, 
if interpreted in the same spirit as above, indicates that rather large relative fluc- 
tuations of temperature, of the order of AT/T ~ 0.12, exist in nuclear collisions. 
It could mean therefore that we are dealing here with some fluctuations existing in 
small parts of the system in respect to the whole system (according to interpretation 
of ||12||) rather than with fluctuations of the event-by-event type in which, for large 



multiplicity N, fluctuations AT/T = 0.06/viV should be negligibly small [10 



We shall now propose a general explanation of the meaning of the function f(x) 
describing fluctuations of some variable x- m paticular, we shall be interested in 
question why, and under what circumstances, it is the gamma distribution that 
describes fluctuations. To this end let us start with the well known equation for the 
variable Xi which in the Langevin formulation has the following form JT^j 



dx 
dt 



- + m 

T 



X = <t> = const > 0. (11) 



Let us concentrate for our purposes on the stochastic process which is defined by 
the white gaussian noise £(£) with ensemble mean 

(£(*)) = (12) 

and correlator (£(£) £(t + At)), which for sufficiently fast changes is equal to 

(£(*)£(* + At)) = 2D5(At). (13) 

Constants r and D define, respectively, the mean time for changes and their variance 
by means of the following conditions: 

(X{t)) = Xo exp (-^j and ( X 2 (t = oo)> = -Dr. (14) 

Thermodynamical equilibrium is assumed here (i.e., t » r, in which case the 
influence of the initial condition xo vanishes and the mean squared of x has value 
corresponding to the state of equilibrium). Making use of the Fokker-Plank equation 

we get for the distribution function the following expression 



ft*) = ~cF~T\ exp 



Mx!) 



rx 

2 / dx 
Jo 



(16) 
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where the constant c is defined by the normalisation condition for f(x)'- Jo°° dxfix) = 
1. K% and K 2 are the intensity coefficients which for the process defined by eq. (|ll| 
are equal to (cf., for example, flTEfl): 





T 

K 2 ( X ) = 2D X 2 . (17) 
It means therefore that as result we have the following distribution function 

-. / \ a— 1 

fix) 

which is nothing but a gamma distribution of variable 1/x depending on two pa- 
rameters: 

(i = — and a = — — . (19) 
D t D 

Returning to the g-notation (cf. eq. @) we have therefore 

q = 1 + rD, (20) 

i.e., the parameter of nonextensitivity is given by the parameter D describing the 
white noise and by the damping constant r. This means then that the relative vari- 
ance lu(1/x) of distribution (pf) is (as in eq. (§)) given by tD. 

As illustration of the genesis of eq. ([□]) used to derive eq. (|20|), we turn once 
more to the fluctuations of temperature |L0|, [II], |12|| discussed before (i.e., to the 
situation when \ — T). Suppose that we have a thermodynamic system, in a small 
(mentally separated) part of which the temperature fluctuates with AT ~ T. Let 
describes stochastic changes of temperature in time. If the mean temperature 
of the system (T) = T then, as result of fluctuations in some small selected region, 
the actual temperature T' equals 

T' = T - b£(t)T, (21) 

where the constant b is defined by the actual definition of the stochastic process 
under consideration, i.e., by £(t), which is assumed to satisfy conditions given by eqs. 
djjD and (|13|). The inevitable exchange of heat between this selected region and the 



rest of the system leads to the equilibration of the temperature. The corresponding 
process of heat conductance is described by the following equation fl7| 



dT 

CpP —-a(T'-T) =0, (22) 
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where c p , p and a are, respectively, the specific heat, density and the coefficient of the 
external conductance. Using T' as defined in (|2~TD we finally get the linear differential 
equation ( |TT| ) for the temperature T with coefficients: r = = -f-T$ = T /r 
and b = t: 



dT 

1H + 



(I 



+ —&£(*) 



T = T . (23) 



c p p c p p 

This result demonstrates clearly that one can think of a deep physical interpretation 
of the parameter q of the corresponding Levy distribution describing the distribu- 
tions of the transverse momenta mentioned before. In this respect our work differs 
from works in which G q (x) is shown to be connected with G q= i(x) = g(x) by the 
so called Hilhorst integral formula (the trace of which is our eq. @) |L3|, [H| but 
without discussing the physical context of the problem. Our original motivation was 
to understand the apparent success of Tsallis statistics (i.e., the situations in which 
q > 1) in the realm of high energy collisions. 

To summarise: if fluctuations of the variable \ can be described in terms of the 
Langevin formulation, their distribution function f(l/x) satisfies the Fokker-Plank 
equation and is therefore given by the Gamma distribution in the variable 1/x- Such 
fluctuations of the parameter 1/x in the exponential formula of physical interest, 
g(x/x), lead immediately to a Levy distribution G g> i(x/x) with q parameter given 
by the relative variance of the fluctuations described by /(1/x)- It should be stressed 
that in this way we address the interpretation of only very limited cases of appli- 
cations of Tsallis statistics. They belong to the category in which the power laws 
physically appear as a consequence of some continuous spectra within appropriate 
integrals. It does not touche, however, a really hard case of applicability of Tsallis 



statistics, namely when zero Lyapunov exponents are involved Hl9[ . Nevertheless, 
this allows us to interpret some nuclear collisions data in terms of fluctuations of 
the inverse temperature, providing thus an important hint to the origin of some 
systematics in the data, understanding of which is crucial in the search for the new 
state of matter: the Quark Gluon Plasma @, [TT[ . 
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